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PREFACE 


The  present  report  is  part  of  a  continuing  effort  to  find 
improved  methods  for  ground-based  detection  of  nuclear  explosions 
conducted  at  high  altitudes  and  in  space. 


SUMMARY 


The  electromagnetic  signal  generated  by  interaction  of 
nuclear— explosion^nduced  electric  currents  with  the  geomagnetic 
field  can  serve  as  a  basis  for  detecting  such  explosions  at  extreme 
distances  from  the  earth.  In  the  present  paper,  an  earlier  theory 
of  this  phenomenon  developed  by  the  authors  is  extended  to  allow 
estimates  of  the  maximum  detection  range.  It  is  concluded  that 


this  electromagnetic  signal  may  be  detectable  at  distances  of  the 
order  of  10  V  Y  kilometers  from  a  nuclear  explosion  of  total 


yield 


Y  kilotons. 

xV  1 
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I.  INTRODUCTION 

In  a  previous  paper  the  authors  considered  the  detection  of 

electromagnetic  radiation  from  a  nuclear  explosion  in  space  caused 

by  Compton— recoil  electrons  from  gamma  rays  scattered  in  the 

materials  of  the  nuclear  device  or  in  a  surrounding  material  shield — 

the  so-called  case  signal.  It  was  concluded  that  the  case  signal 

4  —3/2 

had  a  peak  electric  field  strength  of  about  10  R  v/m  at  a  dis¬ 
tance  R  kilometers  from  the  explosion.  The  maximum  detectable 

0 

range  of  this  signal  was  estimated  to  be  10  km,  independent  of  the 
explosion  yield. 

2 

More  recently  ,  the  authors  considered  another  mechanism 
for  the  generation  of  electromagnetic  radiation  from  nuclear  explo¬ 
sions  which  might  be  useful  for  detection.  This  mechanism  was 
based  on  the  interaction  of  nuclear  explosion— induced  electric  currents 
with  the  geomagnetic  field.  The  theory  of  this  mechanism  was  applied 
to  explosions  which  take  place  well  within  the  earth's  atmosphere. 

It  was  indicated,  however,  that  this  same  geomagnetic  interaction 
mechanism  applies  to  high— altitude  explosions  where  it  can  lead  to 

^W.  J.  Karzas  and  R.  Latter,  Phys.  Rev.  126,  1919  (1962). 

2 

W,  J.  Karzas  and  R.  Latter,  Journ.  of  Geoph.  Res. 

4635  (1962).. 
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quite  intense  electromagnetic  signals  at  great  distances  from  the 
burst  point — thereby  suggesting  its  usefulness  for  detecting  nuclear 
explosions  in  space.  In  the  present  paper,  this  earlier  theory  of  the 
geomagnetic  interaction  mechanism  will  be  extended  so  as  to  permit 
an  estimate  of  the  intensity  of  these  electromagnetic  signals  observed 
on  the  earth's  surface  from  nuclear  explosions  conducted  at  the  great 
distances  of  interest  to  the  detection  problem.  It  is  shown  that  this 
signal  can  be  detected  at  significantly  greater  distances  than  the  case 
signal. 


-- 
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n.  GENERAL  CONSIDERATIONS 


For  nuclear  explosions  below  an  altitude  of  100  or  so  kilo¬ 
meters  the  primary  source  of  the  electromagnetic  signal  is  Compton— 
recoil  electrons  produced  by  explosion  gaunma  rays  interacting  with 

3 

the  atmosphere.  Above  100  or  so  kilometers,  x  rays  from  the  heated 
materials  of  the  nuclear  device  generate  photoelectrons  which  also 
contribute  a  radiating  current. 

The  currents  produced  by  the  gamma  rays  and  x  rays  from 
high— altitude  explosions  have  maxima  in  two  regions.  One  region  is 
near  the  explosion  where  the  gamma— ray  and  x-ray  intensities  are 
highest  even  though  in  this  region  the  air  density  may  be  quite  low. 

The  second  region  is  in  the  earth's  atmosphere  at  altitudes  of  about 
20  to  40  or  so  kilometers  for  gamma  rays  and  of  about  70  to  110  or 
so  kilometers  for  x  rays.  These  altitudes  correspond  to  the  region 
of  peak  absorption  for  explosion  gamma  rays  which  have  a  mean 
energy  of  about  one  mev  and  for  explosion  x  rays  which  have  a  Planck 
distribution  corresponding  to  a  nominal  temperature  of  about  one  kev. 

In  our  treatment  of  the  detection  problem,  radiation  from  the 
currents  near  the  burst  point  will  be  neglected.  For  explosions  at 

3 

At  lower  altitudes,  the  range  of  the  x-ray  photoelectrons  is 
too  short  to  contribute  an  appreciable  current.  Moreover,  the  x  rays 
are  confined  to  a  relatively  small  region  about  the  burst  point,  which 
region  is  made  conducting  by  ionization  from  the  earlier  emission  of 
gamma  rays  and  so  cannot  radiate. 


altitudes  which  are  not  too  great,  the  ionization  in  the  atmosphere 
produced  by  the  gamma  rays  and  x  rays  is  sufficient  to  absorb  this 
radiation  before  it  reaches  the  earth's  surface.  For  explosions  at 
extreme  altitudes,  where  the  gamma^'ay  and  x— ray  ionization  is 
negligible,  this  radiation  can  reach  the  earth's  surface,  but  it  may 
still  be  neglected,  since  its  intensity  turns  out  to  be  substantially 
less  than  that  from  the  atmospheric  currents.  This  lower  intensity 
results  from  the  very  long  mean  free  path  for  production  of  currents 
near  the  burst  point,  which  more  than  compensates  for  the  closeness 
of  this  region  to  the  burst  compared  to  the  great  distance  of  the 
atmosphere.  Moreover,  the  intensity  of  the  signal  from  the  region 
near  the  burst  point  is  further  degraded  by  dispersion  through  having 
to  traverse  the  low  density  ionized  space  medium  for  great  distances 
and  then  having  to  traverse  the  high  density  ionized  ionosphere. 

Figure  1  shows  schematically  the  geometry  of  the  explosion. 
Since  the  currents  contributing  to  the  electromagnetic  field  at  the 
observation  point  are  confined  to  a  relatively  small  region  near  that 
point,  this  geometry  may  be  simplified  by  assuming  that  the  geomag¬ 
netic  field  is  uniform  (See  Fig.  2).  Spherical  polar  coordinates  will 
be  used.  The  origin  will  be  located  at  the  burst  point  and  the  polar 
axis  in  the  direction  of  the  geomagnetic  field  appropriate  to  the 
particular  current  source  being  considered. 
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III.  ELECTRON  CURRENT  AND  DENSITY 

The  number  of  gamma  rays  or  x  rays  emitted  by  a  nuclear 
explosion  per  unit  time  is 

N(t)  -  X  f(t)  ,  (1) 

where  Y  is  the  yield  of  the  explosion  in  the  form  of  gamma  rays  or 
X  rays,  E  is  the  mean  gamma— ray  or  x—ray  energy,  and  f(t)  ex¬ 
presses  the  time  variation  of  the  gamma— ray  or  x-^:ay  emission, 
normalized  so  that 

m 

j  dt  f{t)  «  1  .  (2) 

J 

—  «i 

The  rate  at  which  primary  electrons  (Compton  recoils  and  photoelec— 
tronsj  are  produced  at  a  distance  r  in  direction  e,cp  from  the  explo¬ 
sion  is 


HpriCr,  t)  »  g(i;J  f  (t  -  r/c)  , 


where 


g(L) 


E 


4Tr^ 


(3) 


(4) 


and  the  integral  extends  from  the  burst  point  to  the  point  under  con¬ 
sideration.  X  denotes  the  mean  free  path  for  production  of  primary 


electrons 


“  3  3C  10^  p  /p  cm,  (5) 

V  VO  o  o 

X  *X  plp~^E^  pip  cm,  (6) 

X  xo  '^o  4  X  '^o  ' 


where  P^/p  is  the  ratio  of  sea— level  air  density  to  the  air  density 
at  the  point  of  interest,  and  is  the  x-ray  energy  in  kilovolts. 
The  primary  electrons  are  slowed  down  by  the  atmosphere  and  are 
stopped  in  a  distance  given  by  their  range,  R  ,  where 


2 

R  *R  p/p<s3xl0  p/p  cm, 
Y  ov  '^o  o 

—3  2 

R  •  R  p  /p  «  10  E  p  /p  cm. 

X  QX'^O  X  '^o 


(7) 

(8) 


For  simplicity  we  shall  assume  that  the  electrons  maintain  their 


initial  speed,  v^  ,  throughout  their  range  and  then  abruptly  stop. 

At  a  particular  point  in  space  and  time,  the  primary  electron  density, 
secondary  electron  densities,  and  currents  are  determined  by  adding 
the  contributions  from  all  primary  electrons  made  earlier  and  which 
originated  at  appropriate  positions  and  times.  The  detailed  calcula¬ 
tion  is  given  in  the  Appendix.  We  give  here  only  the  results. 

In  the  absence  of  a  magnetic  field,  the  primary  electron 
density  at  position  r  and  at  time  t  is 


"pri^^  t)  g(r) 


R/v 

r  V  -1 

i  r  °dT'f 

T_T- 

j 

0 

v.  c  y 

(9) 
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where  t  »  t  —  r/c  .  If  we  assume  that  the  secondary  electrons  are 

made  at  a  uniform  rate  during  the  slowing  down  of  the  primaries, 

then  their  rate  of  production  is  obtained  from  the  density  of  primaries 

V  q  E  . 

by  multiplying  by  ,  where  q  =-  and  E  .  is  the  initial 

I-  ^  o  j  j  -I  33ev  pri 

energy  of  a  primary  electron.  Thus 

V  q  R/v  r 

n  (r,t)«-|^  g(r)  j  °dT'f  /c)  ,  (10) 

sec  —  K  —  J  o 

o'-  J 


and  the  secondary  electron  density  is  then 

V  q  T  R/v  r 

n  (r,  t)  Hi g(r)  dr'  I  °dT''f  /c)  ,  (11) 

sec  -  K  J  J  o 

—  a>  O  ^ 


With  no  magnetic  field,  the  net  electron  current  is  in  the  radial 
direction  and  is 

R/v  r 

~  J  T-T'(l-v  /c)  ,  ( 

o  i- 


where  <  equals  the  mean  cosine  of  the  angle  of  scatter  and  is  the 
effective  fraction  of  primaries  moving  in  the  radial  direction.  For 

4 

the  photoelectric  effect,  ‘  ~  T  v  /c  and  for  Compton  scattering 

D  O 

C  W  1  . 

In  the  presence  of  a  magnetic  field,  there  are  two  modifica¬ 
tions  to  the  above  results.  First,  currents  are  generated  in  the 
transverse  as  well  as  the  radial  direction.  Secondly,  as  the  electrons 
are  accelerated  by  the  magnetic  field  the’.r  radial  component  of 


velcx^ity  is  reduced  which,  in  the  gamma-^'ay  case,  causes  the 


current  to  lag  in  time  behind  the  gamma-ray  pulse  and  which,  in  the 
X— ray  case,  causes  an  additional  lag  that  is  a  small  increment  to  the 
lag  already  caused  by  the  low  initial  photoelectron  velocities.  The 
modified  results,  applicable  both  to  the  gamma  rays  and  to  the  x  rays, 
are 


R/v 


1  t  %  a  M  I  o  •  SLn  (jJT 

T  — (1—3  cos  0)t  +psin  0 


(0 


R/v 

j  (r,  t)«-e<v  g(r)  [  ^diT 

O  J 

o 


T  —  (1  — 3  COS  0)  T' 


^  sinwT'  2_  ^ 

+  3  sin  0 -  COS  0  +  Sin  0  cosw 


-e<v^g(r) 


HI 

y  “dT'fj^T- 


(1—3  cos  0)  T* 


^  -  .2-  sinuT' I  .  -  „  .  .  ,x 

+  3  sin  0  — - —  1  sin  0  cos  0  (cosu  t  —  1)  , 


^/v  r  2 

j  (r,t)«  — e<v  g(r)  ax' f  t- (1—3  cos  0)t' 

°  “  o'*  I 

+  3  sin^0  j  sin  0  sinw  t*  , 


(14) 


(15) 


(16) 


where  for  gaunma  rays,  *  »  1  and  3  “  *  rays, 

4 

«  »  —  3  ,  and  terms  in  the  brackets  involving  3»v/c<<l 

5  o 

can  be  negle’cted.  0  is  the  angle  between  the  magnetic  field 
and  the  direction  of  motion  of  the  x  rays  and  gamma  rays, 


—9— 

and  w  is  the  Larmor  fi  equency  for  the  primary  electrons, 
namely, 


eB 


1/2 


w  • 


Y  me 


(l- 

me  V  o  y 


(17) 
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ry.  ELECTROMAGNETIC  FIELDS  FROM 
fflGH-ALTITUDE  CURRENTS 


Equations  (13)  to  (16)  show  that  a  nuclear  explosion  produces 
a  current  and  electron  density  in  the  atmosphere  which  expands  as 
a  pulse  at  the  speed  of  light  radially  from  the  source.  The  magnitude 
of  the  pulse  is  proportional  to  g(r)  of  Eq.  (4)  which  is  large  both 
near  the  source  and  in  the  earth's  atmosphere,  where  most  of  the 
radiation  is  absorbed.  As  remarked  before,  we  shall  limit  our 
considerations  to  the  pulse  produced  in  this  latter  region.  The  time 
behavior  of  the  pulse  is  that  of  a  function  which  varies  rapidly  with 
the  retarded  time  t  •  t  —  r/c  ,  where  r  is  distance  from  the 
explosion.  It  is  convenient,  therefore,  to  rewrite  Maxwell's  equa¬ 
tions  in  terms  of  the  retarded  time. 

In  the  conventional  form.  Maxwell's  equations  are 


7  X  B 


7  X  E 


c  3t 


cat 


7 .  B  «  0  , 


7  .  E  *  4tp  . 


(18) 

(19) 

(20) 
(21) 
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Using  the  transformations  to  retarded  time  t  »  t  —  r/c  ,  namely, 


a  a 
at  ""aT  ’ 


1  a 

-  V  -  e  -  —  , 
-t  — T  — r  c  at 


(22) 

(23) 


where  e^  is  a  unit  vector  directed  radially  from  the  explosion  point, 
Maxwell's  equations  become 


1  1  ^ 

~  1—  .4t. 

—  c— r  at  cat  c-^ 


1 

cat 


1 

-7xE+-^x— ,  , 


1 

7  .  B  -  -  -  0  , 

-  c  at  ’ 


l^^r 


(24) 

(25) 

(26) 

(27) 


As  usual,  j_  represents  the  total  current,  including  "back  currents" 
*or  polarization  currents  in  the  medium.  Conservation  of  charge, 
which  follows  from  these  equations,  becomes 


i£ 
a  T 


+  7  . 


^  cat 


-  0  . 


(28) 


It  is  convenient  to  eliminate  E  from  Maxwell's  equations 
and  get  equations  involving  only  B  ,  and  conversely  .  Doing  this,  we 
obtain  the  following  equations 
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(7^B)  -  ^(vxj)  +11; 

—  r  c  r  c  St 


-  ^  (rB  ) 
r  Sr  r 


0, 


(29) 


Ir,^x3\  4*-  .  v  ^  2  a  f  1  a  .  _  .  2ir 

-  B)q  -  -  (7X^)^  +  _  _  _  (rB^  -  - 


4 


0  ,  (30) 


^  I  "T  f  “ 

—  9  c  -^9  caT|_r 


a  ,  _  .  2t  . 


0,  (31) 


2  2 
where  (v  B)  is  the  radial  component  of  7  B  >  etc.  Analogously, 


—  r 


(7^E)  +  ^  (7  .  j)  +  4t  (7p)  +  7  |- 

—  r  c  r  c  St 


-■|-  (rE  ) 
r  ar  r 


*  0  ,  (32) 


(72E)g^4T(7p)g*||^ 


1  a  /  -  X  ^  2t  . 

F  a?  T  ^0 


0 , 


(33) 


+  4,  (vp)_^  4  I  A 


1  3  /  x:.  \  u.  2t  . 

r  dr  ^  c 


0  ,  (34) 


where 


^  ,  1—L 

St  c  St 


-  j 


(35) 


We  now  observe  that  the  variation  of  the  currents  with  dis¬ 
tance  r  is  much  slower  than  the  variation  of  the  currents  with  c  t. 

This  latter  variation  extends  over  times  t  much  less  than  a  micro— 

4 

second  ,  corresponding  to  the  emission  times  of  the  gamma  rays  and 
X  rays  from  the  explosion.  The  spatial  variation  of  the  currents 
extends  over  distances  like  the  atmospheric  scale  height  ( ~  10  km), 

4 

Latter,  R.  ,  R.  F.  Herbst  and  K.  M.  Watson,  Ann.  Rev.  of 
Nuc.  Sci.  n,  37  1  (1961). 
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which  corresponds  to  times  of  tens  of  microseconds.  For  the  transverse 
currents,  the  resultant  electromagnetic  fields  have  a  rapidly  varying 

5 

time  character  similar  to  the  current  sources  .  Thus,  the  time 
variation  of  the  fields  from  the  transverse  currents  is  much  more 


rapid  than  their  spatial  variation.  We  are  thus  led  to  an  approxima 
tion  in  which  we  keep  only  the  last  bracket  in  Eqs.  (30),  (31),  (33) 
and  (34).  That  is. 


at 


F  h  <^^9)  ^  T  je 


~  0  , 


_a_ 


1^ 

r  a  r  T. 


0  , 


(36) 

(37) 


and  similarly  for  the  magnetic  fields  and  .  For  the  radial 
equations,  to  a  similar  order  of  approximation,  we  obtain  a  non¬ 
radiating  field  given  by 


1  i- 

r  ar 


as 


0  , 


and 


aB  o  aTsin  9j  j 

1  ±_  f  _J£r\  _ 

r  ar  \  at  J  r  sine  ae 


0  , 


(38) 


(39) 


^For  the  radial  currents,  the  radiated  field  does  not  have 
a  time  character  similar  to  the  current  sources.  Rather  its  char¬ 
acter  is  determined  by  the  spatial  extent  of  the  current  source. 

For  the  currents  generated  in  the  atmosphere,  the  spatial  extent 
of  the  currents  is  very  great  and  therefore  the  radial  components 
contribute  a  very  low  frequency  electromagnetic  signal.  This  signal 
is,  in  general,  weak  relative  to  that  due  to  the  transverse  currents 
and  therefore  may  be  neglected. 
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where  we  have  approximately  P~  ~  Jj.  from  Eq.  (28)  .  Integrating 

e 

these  equations,  we  arrive  at  the  final  approximate  equations 


8E 

— ^  +  4Tj  =  0  , 
■'r  ’ 


IS  ,  I.  1  *  21  . 


=  0, 

*  0, 


(40) 

(41) 

(42) 


with  the  corresponding  results  for  the  magnetic  field.  The  first 
equation,  Eq.  (40),  is  the  familiar  one  for  the  radial  electric  field. 

We  shall  not  be  concerned  with  this  non-radiating  component,  but 
only  with  the  propagating  transverse  fields. 

We  recall  that  the  current  density  j  is  to  include  the  second¬ 
ary  currents  in  the  plasma  as  well  as  the  primary  electron  currents. 
Thus, 

i  %ri  *j«ec  ■ 

To  determine  i  we  use  the  Lorentz  model.  We  assume  that  the 

plasma  electrons  are  created  at  rest  and  that  their  motion  is  determined 

only  by  the  local  electric  field  and  by  their  collisions  with  the  air. 

Thus,  if  x(t)  is  the  displacement  of  a  plasma  electron  from  its  initial 

position,  which  we  assume  small  compared  to  important  wavelengths, 

x+v^x--|-  E.  (44) 

—  c  —  m  — 

0 

Results  analogous  to  Eqs.  (41)  and  (42),  but  in  a  planar  approxi¬ 
mation  where  the  factors  of  r  are  not  included,  have  previously  been 
obtained  by  C.  Longmire  (unpublished  1964). 
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For  an  electron  created  at  time  t^  ,  the  current,  —  ^  ,  is  given  by 


2  t  -v(t-t') 

—  ex(t;t  )  »  —  r  dt'e  ^  E(t')  .  (45) 

—  o  m  — 

o 


If  n  (t  )  is  the  rate  at  which  the  electrons  are  made,  then  the 
sec  o 

current  density  is  given  by 


t 

i(t)  »  -  r  dt  n  (t  )  ex(t;t  )  (46) 

^  '  Jo  sec  0  —  0 

m 


or,  from  Eq.  (45)  , 


2  t  -v(t-t') 

i(t).S-Jdfe  "  (47) 

—  OD 


Clearly,  if  the  plasma  electrons  were  disappearing,  as  well  as  being 
created,  the  same  final  expression  would  result — ''sec^^^  then  being 
the  net  electron  density.  Thus,  we  are  led  to  the  following  result 
for  the  plasma  current 


j  ,  (r, t) 

-■'plasma  — ’ 


2 

e 

m 


t  -v(r)Ct-t'] 

J  dt'e  ng^^(r,t') 


(48) 


If  we  substitute  this  expression  for  the  plasma  current  into 

the  field  equations,  we  obtain  for  and  E . 

0  <p 


-1&- 


_b_ 

hr 


rE^(r,T) 


2r  -v(r)(T-T') 


(49) 


hr 


2r  -'J,(£)(t — ’) 


,  2t  .  .  .  „ 

-^T  rj  (r,  t)  »  0 
c 


(50) 


and  similar  results  for  B  and  B  .  These  equations  are  not,  in 


e 


<p 


general,  soluble  because  of  the  spatial  dependence  of  (r)  in  the 
exponential. 

However,  if  v^(r)  is  large  compared  to  the  frequencies  in 


E  and  n  ,  the  equations  simplify  to 

S  60 

3  ,  2x  _  2x  .  „ 

a?  (rE)-f  —  crE  — rj.O, 


(51) 


where  it  is  understood  that  E  -  E  and  i  “*  j  „  •  This  equation 

9,  cp  '  9,  (p 

is  readily  soluble  and  yields 


2,  1  J  -  fdr"G(r”,T) 

E(r,  t)  *  ~  ~  ~  J  dr'r' j(r',T)e  r*^ 
o 


(52) 


where  a(r,t)  is  given  by 
2 


a(r,  t) 


mv  sec  — 
c 


(53) 


In  the  opposite  extreme,  for  a  collisionless  plasma  with  '^(r)  =  0  , 
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the  equations  simplify  to 

2,e^  ^ 


^  (rE)  +  r  [  dT'n  (r,  t')  E(r,T')  +—  r  j(r,  t) 

3r  me  J  sec  —  —  c  — 


0,  (54) 


which  has  the  solution 


E(r,  t)  ^  J  J  ^  (£.'>T') 


o  —  ® 


8x6 

me 


2  T 


1/2. 


r' 


(55) 


Analogous  results  for  the  magnetic  field  are  readily  derived.  In 

particular,  it  is  easily  shown  that  B  »  E.  and  B.  *  -  E 

(p  B  9  Cp 

Thus,  to  obtain  the  transverse  electric  field  from  a  pulse  of 
transverse  current  moving  with  light  speed,  we  add  all  the  contri¬ 
butions  along  the  line  from  the  explosion  to  the  observer.  In  the 
conducting  plasma  these  contributions  are  reduced  by  absorption 
between  the  source  and  the  observer;  in  the  collisionless  plasma, 
dispersion  is  introduced  by  the  plasma. 


If  is  independent  of  position,  it  is  easy  to  show  that  the 

solution  of  Eqs.  (49)  and  (50)  is  obtained  from  Eq.  (55)  by  replacing 

'"c  ^  ^cl 

E(r,  t)  with  E(r,  t)e  and  j(r,  t)  by  j(r,  t)e 
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V.  DETECTION  CONSIDERATIONS 

It  is  clear  from  Eqs.  (52)  and  (55)  that  the  electromagnetic 
field  observed  from  high— altitude  nuclear  explosions  depends  only 
upon  the  distribution  of  current  and  electron  density  along  the  line 
between  the  burst  point  and  the  observation  point.  However,  because 
of  the  exponential  variation  of  atmospheric  density,  the  magnitude  of 
the  electromagnetic  field  depends  quite  sensitivity  upon  just  where 
the  currents  and  electron  densities  are  generated  in  the  atmosphere. 
Moreover,  the  occurrence  of  the  conductivity  a  (£,  t)  in  the  exponen¬ 
tial  of  Eq.  (52)  indicates  that  the  magnitude  of  the  electron  density 
at  each  position  in  the  atmosphere  is  also  quite  important.  The 
result  is  that  the  magnitude  of  the  observed  electromagnetic  field 
is  quite  sensitive  to  the  gamma^ray  and  x— ray  yield  of  the  explosion, 
the  distance  between  the  observation  point  and  the  burst  point,  the 

temperature  of  the  explosion,  and  the  angle  of  incidence  of  the 

0 

radiation  on  the  atmosphere  .  A  comprehensive  description  of  the 
electromagnetic  field  would  therefore  require  a  parametric  study 
for  beyond  the  intent  of  this  paper  which  is  simply  to  indicate  the 
possible  usefulness  of  this  signal  for  detection.  What  we  shall  do 

0 

Effects  due  to  the  earth's  ionosphere  are  not  important  in 
the  present  considerations,  since  the  currents  which  generate  the 
electromagnetic  signal  are  produced  below  the  ionosphere.  For  a 
detailed  discussion  of  the  deposition  of  explosion  x  rays  and  gamma 
rays  in  the  atmosphere,  see  R.  Latter  and  R.  E.  LeLevier,  Journ. 
Geophys.  Res.  68,  1643  (1963). 


here  is  simply  indicate  the  crude  order  of  magnitude  of  the  electro¬ 
magnetic  field  strength  under  typical  explosion  conditions  and  then 
indicate  roughly  the  limiting  detection  distance. 

First,  we  shall  consider  explosions  which  are  not  too  distant 
from  the  earth,  which  have  a  temperature  of  about  one  kev,  and 
for  which  radiation  is  incident  almost  normally  on  the  atmosphere, 
^ecifically,  these  conditions  mean  that  the  ionization  produced 

deep  in  the  atmosphere  by  the  x  rays  and  gamma  rays  is  high,  so 

g 

that  in  the  regions  of  interest  is  large,  allowing  the  approxi¬ 
mation  Eq.  (52)  to  be  applied,  and,  moreover, 


c 


r 

j  dr'  a  (r*,  T  )  >  1  . 
o 


(56) 


This  condition  implies  for  x  rays  that  a  one  kiloton  explosion 

4 

IS  within  about  a  few  times  10  km  from  the  earth's  surface 

2 

and  for  gamma  rays,  within  about  a  few  times  10  km. 


For  low  pressures,  large  electric  fields  and  electron-volt 
energies,  the  collision  frequency  of  electrons  in  air  is  approximately 
independent  of  the  E— field  and  has  a  value  v  ■  4  x  10^^  p/p  sec“^  ; 
see  R.  A.  Nielsen  and  N.  E.  Bradbury,  Phys.  Rev.  51,  69°(1937). 
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To  estimate  the  magnitude  of  this  radiation,  we  first  simplify 
Eq.  (52).  We  note  that  the  factor  r'  in  the  integrand  of  Eq.  (52)  can  t 
replaced  by  r,  the  distance  to  the  observation  point,  since  the 
important  currents  ar^  confined  to  a  thin  layer  in  the  atmosphere 
near  the  observation  point.  Thus 


dr'j  (r',T)e 

e,  {p 


r' 


dr"a(r",  T) 


(57) 


Now  we  note  that  for  r'  sufficiently  near  r  , 


dr"a(r",T)«  0  and  j(r',T)ss  0,  (58) 

since  the  gamma  rays  and  x  rays  are  rapidly  absorbed  by  the  atmos¬ 
phere  below  some  low  altitude.  As  r*  decreases  and  the  altitude 
increases,  we  finally  reach  a  radius  r  at  which  the  x-ray  and 
gamma-ray  ciirrents  and  ionization  increase  rapidly.  Thus  we  can 
rewrite  Eq.  (58)  as 


E  (r,  T)  2: 

«,  cp 


2r 


r 

r  o 


J 

o 


dr' 


jfl  (r’,T)e 

8,  q) 


_  ^ 
c  J. 


dr"a(r",T) 


(59) 


or  making  an  obvious  change  of  variables, 


dr"a(r",  t) 


L  (r'(x),  T) 

d  X  -  ^  -*8.  ■ 

CT(r'(»),  t) 


-X 


(60) 


jg  (r'(x),T)  _  I 

It  is  easy  to  show  that  — varies  slowly  compared  to  e  * 


and  therefore 


E.  (r,T) 

a,  CP 


a(r  ,  T)  '0  ®  °  }  ’ 

o 


jfi 

0,  cp  o 


a  (r^ ,  T) 


Substituting  Eos.  (15),  (16)  and  (53)  into  Eq.  (61),  we  find 
« mu  V  R  (t) 

.  _ E_  Qjg  ■■  (62) 

P  eq  N(T)  ’  ^ 


where 


T  R  /v  r 

■KT/  ,  ^  o  .  / ,  a  2„.  ,  ,  -  .  sinu  T' 

N(t)  *uj  dT  j  dT'f  T  —  (1  — 6cos  0)T'+psin  9  — - —  , 

-*  o  L 

J„(t)  »  1^  °dT' f  T  —  (1  — 6  cos^9)t' +9sin^9  8in9cos9  (c 

0  J 

o  L  .  J 


V"'  ■  r 


sin^9  sin9cos9  (cosut'— 1)  ,  (64) 


°dT'f  T— (1— 0  cos^9)t' +3sin^9^^^^^^  sin9sinuT'.  (65) 

b> 


We  note  that  for  gamma  rays,  which  deposit  most  of  their  energy  in 
the  altitude  range  from  20  to  40  km. 


R  *“7  0 

—  ~  2  X  lU  to  3  X  10  sec. 

V 

o 


Equation  (61)  gives  the  result  obtained  by  a  more  direct 
calculation  in  the  paper  cited  in  footnote  2. 
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4  -7 

Since  the  gamma  rays  are  emitted  in  times  of  about  10  sec,  we  see 

that  order— of— magnitude--wise  f(T)  can  be  replaced  by  a  delta  function. 

4 

Similarly,  for  x  rays  which  are  emitted  in  times  considerably  less  than 

g 

10  sec,  and  are  absorbed  in  the  altitude  range  from  70  to  110  km, 

— ^  ~  10  to  5  X  10  ^  sec,  (67) 

o 

_2 

where  R  ~  10  electrons  with  energies  of  a  few  kilovolts, 

so  that  again  f(T)  can  be  replaced  approximately  by  a  delta  function. 
Thus, 


Jr) 


J Jr)  mt 


X  , 

X 

^  wR 

3c 

wR 

^  mR 

3c  ’ 

^  W~ 

sin  9  cos  0(1-  cos  X) 

- 

-  uR 

2  2 

1  — 3cos  0— 3sin  0cos  X 

,  X 

0  , 

X 

^  wR 

Pc 

sin  0  sin  X 

Y 

.  uR 

2  2 

1—3  cos  0—3  sin  9  cos  X 

9  A. 

3  c 

0  , 

X 

^  «R 

3  c 

where  X  is  determined  from 


2  2 

(1  —  3  cos  0)  X  —  3  sin  9  sin  X  «  w  t  . 


(68) 


(69) 


(70) 


(7  1) 


Numerical  evaluation  of  using  3«v^/c»0.  95  correspond¬ 

ing  to  1  mev  Compton— recoil  electrons  and  3  *  v  /c  »  0.  05  to  0.  1 

o 

corresponding  to  1  to  a  few  kev  x-ray  photoelectrons,  gives  the  peak 
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J  /N~  1/(1— P)  •  Hence,  ordei>-of-magnitude-vrlse 

e,cp 


E 


«  m  w  V  R 
c 


9,cp  eq(l-0) 


(72) 


.-28 


For  gamma  rays  wherein  c  JS  1,  m  *^9  x  10  gm, 
w  «  1.  67  X  10®  cps,  e  -  4.  8  X  10“^°,  q  «  3  x  10^,  and  R  «  3  x  10^ 


pi  p  cm,  we  have 
o 


E  ~  2  esu  or  6  X  10  v/m  . 
0,CP 


(73) 


6  2 

And  for  x  rays,  wherein  t  ~  0.  1,  w  ®  5.  3  x  10  cps,  q  ~  10  , 

_2 

and  R  «  10  p/  , 


E  ~  3.  5  X  10~*  esu  or  10  v/m  . 

0,<P 

Both  of  these  fields  are  sufficiently  intense  so  as  to  be  easily 
detectable. 


(74) 


We  now  turn  to  the  problem  of  estimating  the  intensity  of  the 
electromagnetic  signal  from  explosions  at  extreme  distances  from 


mean  photoelectron  energy  of  about  3  kev  is  assumed, 
since  the  x  rays  which  contribute  most  to  Eg  are  those  which 
penetrate  deepest.  For  the  range  of  low— enei^gy  electrons,  see 
Eq.  (8)  and  L.  Katz  and  A.  S.  Penfold,  Rev.  Mod.  Phys.  28  (1952). 

Actually,  for  the  assumed  close  distance  of  the  explosions 
and  the  assumed  normal  incidence  of  the  radiations  on  the  atmosphere, 
the  pre-ionization  produced  by  the  gamma  rays  is  sufficiently  intense 
to  absorb  the  x-ray  produced  signal  before  it  reaches  the  earth's 
surface.  Thus,  under  these  circumstances,  only  the  gamma— ray 
produced  signal  is  useful  for  detection. 
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the  earth,  and  thereby  obtain  the  limit  of  detection.  For  suffi¬ 
ciently  great  distances,  we  may  neglect  n  and  a  in  Eqs. 

S6  C 

(52)  and  (55),  so  that 


0,cp 


2£ 

c 


(75) 


and  the  integration  is  over  only  the  regions  in  the  atmosphere  near 

the  earth's  sxirface  where  gamma  rays  and  x  rays  are  deposited. 

Substituting  Eqs.  (4),  (15)  and. (16)  into  Eq.  (75)  and  making 

the  delta-function  approximation  for  f(T)  ,  we  find 

r' 


E 


(76) 


where  J.  (t)  is  defined  by  Eqs.  (69)  and  (70).  Since  the  integra— 
e,cp 

tions  extend  over  a  region  small  compared  to  r'  and  since  the 
r  dr ' 

integral  J  over  the  important  region  is  large. 


2E  J.2 


(77) 


Finally,  it  is  clear  that  J.  (t)  ~  1,  so  that 

0>«P 


e,(p 


epc  X 
2E  2  ’ 


(78) 
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where  T  is  the  effective  temperature  in  degrees  Kelvin  of  the 

cosmic  noise  at  the  wavelength  X  in  meters  and  a  is  the  detector 

bandwidth  in  cps.  Since  the  signal  to  be  detected  has  frequencies  from 

one  to  ten  megacycles  (See  Eqs.  (66),  (67),  (69),  (70)  and  (71)^,  the 

7  2 

detector  bandwidth  should  be  adjusted  to  a  ~  10  cps  at  X  10  m . 
Under  these  circumstances  the  effective  temperature  of  the  cosmic 
noise  background  is  T  ~  10  °K.  Hence,  the  interfering  noise  power 
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is 

P--  ^  1.  7  X  10  w/m^  .  (82) 

N 

For  reliable  detection — using  coincidence  techniques — we  require 

the  peak  signal  power  P  to  exceed  the  noise  power  by  a  factor  of 

s 

ten,  that  is, 

P  »  ^  -  10  P„  .  (83) 

s  4t  s  N 

Or 

E  «  8  X  10  esu  or  2. 4  x  10  ^  v/m  .  (84) 

s 

Thus,  from  Eqs.  (79)  and  (80),  the  limiting  detection  range  for 
gamma  rays  is 

and  for  x  rays  is 

■■km  ~ 
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VI.  DISCUSSION 

An  interesting  consequence  of  Eqs.  (72)  to  (74)  and  Eq.  (78) 
to  (80)  is  that  for  near  explosions,  the  gamma— ray  generated  electro¬ 
magnetic  signal  is  more  intense  than  the  x-ray  generated  signal, 
while  for  distant  explosions  the  opposite  is  true.  To  understand  this 
result  we  note  that  for  distant  explosions  the  intensity  of  the  radiated 
field  is  proportional  to  the  yield  of  the  explosion  appearing  as  gamma 
rays  or  x  rays — thereby,  as  would  be  expected,  the  x-ray  generated 
signal  should  be  larger.  As  the  burst  point  comes  nearer  to  the 
earth's  atmosphere,  the  intensity  of  the  signal  begins  to  level  off 
and  becomes  independent  of  yield.  For  xrays,  this  leveling  off  occurs 
at  a  very  much  greater  distance  than  for  gamma  rays — thus  the 
smaller  intensity  of  the  x-ray  generated  signal. 

A  final  point — ^in  concluding  that  the  detection  range  of  the 
x-ray  generated  signal  was  given  by  Eq.  (86)  and  exceeded  the  gamma- 
ray  generated  signal,  it  was  assumed  that  the  x— ray  yield  of  the 
explosion  was  about  one-half  the  total  yield.  Actually,  however, 

2 

this  yield  might  be  considerably  less — perhaps  by  a  factor  of  10  to 
3  4 

10  ,  if  shielding  were  employed  .  In  this  case  the  limiting  detection 

range  would  be  determined  by  the  gamma— ray  generated  signal. 
Shielding  of  the  gamma-ray  signal  might  also  be  possible,  thereby 
degrading  the  detection  range  even  further. 
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APPENDIX 


In  this  appendix  we  derive  the  formulae  for  the  primary  elec¬ 
tron  current  and  secondary  electron  density  as  functions  of  space  and 
time. 

The  rate  of  production  of  primary  electrons  at  a  point  r',  at 
time  t',  with  initial  motion  in  direction  Q  ,  is  given  by 

r' 

)  *  g(r')f(t  -- ^)  h  (Q)  ,  (A-1) 

where  the  angular  distribution  h^)  is  so  normalized  that 

Jdn  h(n)  -  1  ,  (A-2) 

Jdn  h(n)  -  «  ,  (A-3) 

and  is  the  cosine  of  the  angle  between  the  direction  of  motion  of 
the  electrons  and  the  initial  gamma-ray  or  x— ray  direction.  For 
gamma— ray  produced  Compton  electrons,  <  w  1.  For  photoelectrons 

4 

from  X  rays,  <  »  —  v  /c  ~  o.  05  to  0.  1  for  1  to  a  few  kev  electrons. 

0  o 

From  symmetry,  the  other  components  of  ^  integrate  to  zero. 

The  primary  electron  current  and  density  at  a  position  r 
and  time  t  are  obtained  by  adding  the  contributions  from  those 
electrons  made  at  earlier  times  and  appropriate  positions  such  that  their 
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motion  brings  them  to  r  at  time  t  .  That  is, 


j(r,  t)  »  —  e  Idr'  fdt'  dflg(r')f(t' - )h(f2)  v(t,  t';  n) 

J  —  J  J  “  “  C  ““  ““ 


r  “  r  *  “ 


t  n 

'  -  Jdt"  v(t",t';n)  , 
t'  J 


(A-1) 


where  v  (t,t';^)  is  the  velocity  at  time  t  of  an  electron  made  at  t' 
with  initial  direction  Q  .  The  density  of  primary  electrons  is  given 
by  the  same  expression  without  the  factor  —  ev.  In  the  absence  of  a 
magnetic  field  or  any  slowing  down  mechanism,  the  velocity  of  the 
electrons  is  constant,  so  that 


v(t,t';£)  -  v^n  ,  (A-5) 

t 

[  dt"  v(t",  t';n)  •  V  n(t-t').  (A-6) 

t’ 

Since  the  electrons  have  a  finite  lifetime  given  by  their  range  R  divided 
by  their  initial  speed,  the  integral  over  t'  in  Eq.  (A— 4)  extends  only 
over  the  finite  range  from  t  —  R/v^  to  t  .  We  shall  assume  that 
the  distances  the  electrons  travel  are  small  compared  to  r,  which  is 
the  case  for  currents  generated  in  the  atmosphere  by  high— altitude 
nuclear  explosions.  Thus,  the  direction  from  the  source  to  r'  can  be 
taken  to  be  the  same  as  that  to  r  and  the  current  simplifies  to 


i  (r,  t)  »  —  ev  g(r  )  dr'  f  dt'  f(t'  —  —  )  f  dll  n  h(n  ) 

^  °  ^  Ur/v  - - 

o 


r  —  r'— V  (t  —  t')l^ 
o  r 


(A-7) 
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With  t  —  t'  E  T',  t  —  r/c  E  T,  we  have 


R/v  ^  V  t' 

j(r,t)*  — ev  g(r)  [  ° dr*  |  dfi  h(Q)  n  f  fr— T'  +  -^n^.  (A— 8) 

-  o—  J  4  —  —  —  V  cry 

o 


For  Compton  electrons  h(f2)  is  essentially  a  delta  function  in 
the  initial  direction.  Thus  the  current  is  radial  and 


R/v 


g(£)  1  dT'f 

O  J 


R/v 

n  i(r,  t)a  g<£)  J 

^  o 


T-  T'(l - ) 

c 


T-  T'  (1  -  —  ) 
C 


(A-9) 

(A-10) 


For  photoelectrons,  we  expand  the  function  f  in  powers  of 


V  T' 

—  Q  <  —  Q  .  Thus 
c  r  c  r 


V  T'  V  t' 

f(T-T'  +•—  n  )W  f(T-T')  +  -2-  n  f(T-T')  + 

c  r  c  r 


(A-11) 


The  integration  over  n  in  (A— 8)  converts  Eq.  (A— 11)  into  an  expansion 

V  T'«  ^ 

o  iRc 

in  terms  of  -  which  is  bounded  by  —  and  which  is  small  compared 

c  c 

to  one  at  all  altitudes  of  importance.  Thus  for  photoelectrons, 

:/v 


R/v 

i^(r,t)w-etv  g(r)  |  °dT'f(T-T')  , 

r  o  V 


R/v 


^  °dT'f(T-T’). 

pri —  —  j 


(A-12) 


(A-13) 


We  note  that  although  different  approximations  have  been  used. 


we  can  write  the  results  in  one  form  for  both  Compton  electrons  and 
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where  v(t')  is  the  velocity  of  an  electron  at  time  t'  after  it  is  made 
and  x(t')  is  its  displacement  in  the  radial  direction  (direction  of  the 
initial  gamma  ray  or  x  ray)  from  its  initial  position.  As  before, 
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«  w  1  for  Compton  alectrons,  and  we  drop  x(t)/c  in  the  argument  of 
f  for  photoelectrons. 

If  only  a  static  magnetic  field  is  present,  the  motion  of  the 
recoil  electrons  is  very  simple.  The  component  of  velocity  along  the 
magnetic  field  is  unchanged  and  the  motion  in  the  plane  perpendicular 
to  the  field  is  circular  with  angular  frequency 
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t  (1-8  ) 
In  terms  of  spherical  coordinates, 
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where  9  is  the  angle  between  the  radial  direction  and  the  mag¬ 
netic  field.  When  these  expressions  are  substituted  into  Eqs.  (A— 16) 
and  (A— 17),  we  obtain  the  results  given  in  Section  IH,  Eq.  (13)  to  (16). 

We  have  neglected  several  effects  in  deriving  these  expressions. 
First,  we  have  assumed  that  the  electrons  have  a  constant  speed,  and 
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then  abruptly  stop  after  an  appropriate  time,  R/v  Actually,  of 

o 

course,  they  slow  down  as  they  lose  energy  by  ionizing  the  air.  Since 
the  rate  of  energy  loss,  F(v)  ,  is  known  for  electrons,  it  can,  of  course, 
be  included.  We  then  have  to  solve  simultaneously  the  equations  of  motion 
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To  do  this,  we  note  that  a  straightforward  integration  of  Eq.  (A— 24) 


gives 

Y 


p  d  Y  ^ _ ^ 

J  F(y)  * 


(A-26) 


me 


where  y^  •  v(t  »  0)  .  Solving  this  equation  for  v  as  a  function  of  T  we 


obtain  the  electron  speed 
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Further,  letting  v*v^  where  u  is  the  unit  vector  in  the  direction 


of  the  velocity,  we  find 
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Thus  the  electron  velocities  are  the  same  as  found  before,  with  v 


replaced  by  v(t)  •  c 


1  - 


1 


V  ^(t) 


1/2 


and  b>  T  replaced  by  9  (t)  . 
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The  net  effect  is  to  reduce  the  velocity  to  zero  during  the  slowing  down 
and  to  increase  the  Larmor  frequency  u  as  v  decreases  to  one. 

Of  greater  concern  is  the  question — what  magnetic  field  should 
be  used  in  determining  the  electron  motion?  Under  conditions  where 
the  currents  generate  only  weak  fields,  we  need  use  the  ambient  mag¬ 
netic  field  only.  However,  under  conditions  where  the  explosions  are 
not  too  distant,  the  currents  can  give  rise  to  intense  crossed  E  and 
B  fields  which  can  significantly  affect  the  motion  of  the  electrons. 

For  x-ray  photoelectrons,  however,  the  fields  even  under  these  condi¬ 
tions  are  sufficiently  small  so  that  only  the  ambient  magnetic  field  need 
be  used.  For  Compton-^recoil  electrons,  the  fields  are  large  so  that 
these  self-generated  fields  must  be  considered.  To  estimate  the 
effect  of  these  Compton-current  generated  fields,  we  note  first  that 
the  electric  field  produced  by  the  electrons  is  transverse  and  in  such 
a  direction  as  to  oppose  the  transverse  motion  of  the  electrons.  The 
generated  magnetic  field,  however,  is  in  such  a  direction  as  to  add 
to  the  ambient  geomagnetic  field  and  increase  the  transverse  accelera¬ 
tion.  The  intensities  of  the  generated  electric  and  magnetic  fields, 
which  are  equal  in  magnitude,  depend  upon  the  amount  of  transverse 
current  and  are  time  dependent.  It  is  clearly  impractical  to  try  to 
solve  rigorously  the  self-consistent  problem  for  the  electron  motions. 
What  we  can  do  first  is  to  observe  that  the  total  force  on  an  electron  is 

F--e  E  +  -^x(B+B^)j  ,  (A-3 1) 

where  E  and  B  are  the  fields  due  to  the  currents  and  B  is  the 

— o 

ambient  field.  Since  |e1  ■  |b|  ,  the  net  transverse  force  is 
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proportional  to  E  (1  —  — )  plus  the  transverse  component  of  —  x 
If  the  radial  electron  velocity,  v^  ,  is  approximately  light  speed, 
the  effect  of  the  induced  fields  is  greatly  reduced.  For  Compton 
electrons,  then,  the  initial  turning  is  due  primarily  to  the  ambient 
field. 

We  can  next  approximately  solve  the  self-consistent  problem 
for  Compton  electrons  by  taking  a  simple  special  case.  We  let  the 
gamma— ray  source  vary  as  a  delta  function  in  time,  and  take  the 
initial  direction  of  the  electron  motion  at  right  angles  to  the  ambient 
field.  Then  the  equations  of  motion  are 
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where  x  is  the  direction  of  initial  motion  and  y  is  the  transverse 
direction.  For  not-too-distant  explosions,  for  which  the  fields  are 
highest,  we  find  from  Eqs.  (61),  (A-16)  and  (A-17)  that 
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It  is  straightforward  to  show  that  the  maximum  field  strength  occurs 


at  T  »  0  and  has  the  value 
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lf  we  neglect  the  generated  field  in  Eqs.  (A-32)  and  (A-33),  and 
consider  the  motion  only  in  the  ambient  geomagnetic  field,  we  find 
the  peak  value  of  the  electric  field  again  occurs  at  t  ■  0  and  is 
given  by 
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Thus  the  effect  of  the  generated  fields  is  to  decrease  the  peak  field 
by  the  factor 
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where  parameter  values  corresponding  to  1  mev  electrons  have  been 
used. 

This  example  indicates  that  as  long  as  the  peak  radiated  field 
strength  occurs  near  the  beginning  of  the  current  pulse,  the  effect 
of  the  generated  fields  on  the  electron  motion  does  not  significantly 
change  the  order  of  magnitude  of  the  estimate  for  the  peak  radiated 
field.  The  generated  fields  do  affect  the  time  variation  of  the 

radiated  field.  But  since  the  range  of  detection  depends  principally 
on  the  peak  field  strength,  the  effect  of  the  generated  fields  on  the 
electron  motion  may  be  disregaurded. 
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Figure  2 

SIMPLIFIED  GEOMETRY  OF  THE  EXPLOSION 


